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ABSTRACT 

In a recent paper, we found evidence for an increase in the accretion rate during photospheric radius ex¬ 
pansion bursts, quantified by a variable normalization factor f a on the pre-burst persistent emission. Here we 
follow this result up on a much larger sample of 1759 type I X-ray bursts from 56 sources. We show that the 
variable persistent flux method provides improvements in the quality of spectral fits for type I bursts, whether or 
not they reach the Eddington luminosity. The new approach has an estimated Bayes factor of 64 improvement 
over the standard method, and we recommend the procedure be adopted as standard for analysing type I bursts. 

We show evidence that the remaining discrepancies to a formally consistent spectral model are due to the burst 
component deviating significantly from a blackbody, rather than variations in the spectral shape of the persis¬ 
tent emission component. In bursts that do not show radius expansion, the persistent emission enhancement 
does not exceed 37% of the Eddington flux. We use this observation to constrain the Eddington flux of sources 
for which /> ( id has not been directly measured. 


1. INTRODUCTION 

Typ e I X-ray bu rsts, discovered in the 1970s (iGrindlav et all 
1 976: Belian et al.ll!9761) . arise from thermonuclear runaways 
in ac creted H/He material on the surface of a neutron sta r 
(e.g., IWooslev & Taamlfl976t IStrohmaver & Bi Idstenl 120061) . 
Gas flowing from a low-mass stellar companion accum ulates 
on the surface o f the neutron star (e.g., IWooslev & Taaml 
H97(l 2 ossl 119771) . There it undergoes hydrostatic heat¬ 
ing and compression, triggering unstable nuclear burning 
when the temperature and pressure are high enou gh (e.g., 
iFuiimoto et all IT98H IStrohmaver & Bildstenl 12006 ). These 
events are observed as a sudden increase in the neutron star’s 
X-ray lu minosity, to many times the persistent lev el (see re¬ 
views by iLewin & Jossl 119811 ILewin et~akl 1 1 9931) . Typical 
bursts exhibit rise times of 1—10 du ratio ns of a few tens 
of seconds to a few minutes (e.g.. lGalloway et aDl2008 ). and 
emit a total energy of 10 39 —10 40 erg. 

Type I bursts provide a probe into the conditions on the 
surface of the neutron star and yield in sights into its inter¬ 
nal structure (lLattimer & Prakashl 12007 ). Performing such 
measurements requires an accurate determination of the flux 
throughout the burst. This in turn necessitates a good under¬ 
standing of the X-ray spectra since bolometric flux_cannot 
be m easured directly from photon count dBlissett & Cruisel 
Il979l) . and we must also know how to separate the flux from 
the expanding atmosphere from that arising from other sites 
in the neutron star system. One such contribution is the ac¬ 
cretion disc itself. It emits X-rays through the conversion of 
gravitational potential energy of the inspiralling gas. This 
emission is present before the burst, and has traditionally 
been assumed to remain constant (both in spectral shape and 


in intensity) throughout the burst (e.g.. 

van Paradiis & Lewin 

1986 

; Lewin et al. 1993c Kuulkers et al. 

20(E| Galloway et al. 

2008 

). 


iMuno et akl (120001) and IStrohmaver & Brownl (120021) al¬ 
lowed for the persistent emission to be suppressed entirely. 
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by subtracting the instrumental background only, but did not 
find that the spectral fits were improved by doing so. In a 
study of neutron star radii derived from cooling tail spectra, 
iGiiver et~akl (120121 ) noted the possibility of variations in the 
persistent emission spectrum and excluded bursts whose pre¬ 
burst persistent flux was more than 10% of the Eddington flux 
of that source, so that the contr ibuti on from persistent emis¬ 
sion is always small, lin’t Zand et all (120131) studied a PRE 
burst from SAX J1808.4-3658 using combined Chandra and 
RXTE data and found that an observed excess of photons at 
both low and high energies can be well described by allow¬ 
ing a 20-fold increase of the pre-burst persistent emission. 
This approach was app lied to 332 PRE bursts observed with 
RXTE by IWornel et al.l (2013; hereafter WGP13), where we 
found that allowing the persistent emission to vary in intensity 
during a burst improves the spectral fits, and that the persis¬ 
tent emission usually increases temporarily to several times 
its pre-burst level. 

Since the intensity of the persistent e missio n is expected 
to increase with increasing accretion rate. [WGP13I interpreted 
this effect as a temporary accretion enhancement, possibly the 
result of the disk losing angular momentu m via radiation drag. 
Other interpretations are possible; iin’t Zand et al.l 120131) ar¬ 
gue for reprocessing of the burst spectrum in the accretion 
disc. 

The "variable persistent flux" approach has_since been ap¬ 
plied successfully by other studies (e.g. . lin’t Zand et aT1l20 1 3h 
IKeek et al.EoTllPeille et al.l2014l) . and it is becoming clearer 
that bursts affect the accretion disc and can be used to probe 
accretion physics. In an inv estig ation of a long duration burst 
from 4U 1636-536. IKeek et al.l (120141) found that the persis¬ 
tent emission was enhanced by a factor of « 2 and also in¬ 
ferred a significant modification of the inner accretion disc. 
A depletion of the accretion disc through radiation drag dur¬ 
ing a burst was predicted by IWalkerl 1199 2). who found in 
numerical simulations that the accretion rate after the burst 
can decline temporarily to half the pre-burst level. iPeille et afl 
(120141) performed a study of quasi-periodic oscillations from 
4U 1636-536 and 4U 1608-522, and found that these oscil¬ 
lations were suppressed for several tens of seconds after the 
burst. They attribute this phenomenon to a depletion of the in- 
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ner accretion disc, though its fast recovery (5-10 times faster 
than the viscous time) remains a mystery. Clearly a greater 
observational and theoretical understanding of the behaviour 
of accretion onto type I bursters is required. 

Disentangling the burst component of the spectrum from 
the persistent component is difficult, as the various compo¬ 
nents are usually spectrally degenerate. Furthermore, we 
do not have spatial resolution and therefore cannot attribute 
different spectral components to different sites in the neu¬ 
tron star system. This problem will be further compli¬ 
cated if, as we expect, both the persistent and burst com¬ 
ponents can vary in spectral shape during a burst. In¬ 
deed, it has long been known that the spectral shape of 
the persistent emission, as w ell as its intensi ty, is a func¬ 
tion of accretion rate (lHasinger & van der Klislfl989t see also 
iLvu et all 12014 ). There are, however, suggestions that be¬ 
tween bursts the accretion spectrum remains constant in shape 
on timescales longer t han those of a burst, t hough its in¬ 
tensity may vary (e.g.. iThompson e t ah 2005t iBagnoli et al.l 
1201 3b iLinares et al.l 120141) . During a burst it is, of course, 
much more difficult to detect changes in the accretion spec¬ 
tral shape. A sudden loss of high energy (> 30 keV) pho¬ 
tons has been detected from several sources, but this effect 
has been attributed to ra pid cooling of the corona and is not 
related to accretion rat e dMaccarone & Coppil 120031 : iJi et al.l 
120 1 3b IChen et all 120131) . Changes in the shape and intensity 
of the persistent emission during a very long_(20,000s) burst 
from 4U 1636-536 were reported bv iKeek et all ( 2014 ). but 
variations on these timescales can happen in the absence of 
bursts, and the bursts we consider in this paper are very much 
shorter. 

PRE bursts change the structure of the neutron star photo¬ 
sphere, and this is likely to cause the spectrum of the photo¬ 
sphere to deviate from its usual near-blackbody shape. Such 
chang es ar e unrelated to the accretion spectrum. As pointed 
out bv lWGPl 3l these will confound spectral analyses of per¬ 
sistent emission during a burst. Such complications can be re¬ 
moved by considering non-PRE bursts, but care must be taken 
to account for lower sensitivity due to their intrinsically lower 
flux. 

In this paper we investigate the relationship between accre¬ 
tion and type I bursts. We apply the fitting method of WGP13 ; 
to every type I burst detected by RXTE. In addition, we in¬ 
vestigate the variability of the accretion spectra immediately 
before and after bursts, to determine whether our spectral fit¬ 
ting procedure can be confounded by rapid changes in the ac¬ 
cretion spectrum. In ^5] we develop a method for obtaining a 
lower bound on the Eddington flux for sources that have not 
yet been observed to undergo PRE bursts. 

2. DATA ANALYSIS 

We used observational data from the Rossi X-Ray Tim¬ 
ing Explorer (RXTE), publicly available through the High- 
Energy Astrophysics Science Archive Research Centre 
(HEASARCfl The observations date from shortly after the 
satellite’s launch on December 30, 1995 to the end of the 
RXTE mission on January 3, 2012. Ou r samp le of type I bursts 
is based on the burst catalogue of ( iGallowav et al.l[2008l : here¬ 
after G08), with the addition of 880 type I bursts detected 
after the publication of that paper. The entire RXTE burst cat¬ 
alog forms part of the Multi-INstrument Burst Archive (MIN- 

3 See http://heasarc.gsfc.nasa.gov 


BARfJ 

We classified type I bursts as either belonging to the pho- 
tospheric radius expansion (PRE) class of bursts, or not, ac¬ 
cording to the criteria for radius expansion described in §2.3 
of IG08L In brief, these criteria define radius expansion as 
an increase of the surface area of the photosphere together 
with a decrease in its temperature. A small number of bursts 
fulfil some, but not all, of the criteria and are classified as 
"marginal" PRE bursts. 

Our data reduction procedures are the same as those in §2 
of lWGP13l unless otherwise stated, and we refer the reader 
to that work for further details. To estimate and remove the 
instrumental background we used the full-mission, “bright” 
source (>40 counts s -1 ) models released 2006 August 6 with 
the pcabackest tool. Table |T] lists interstellar absorption 
column densities for all burst sources and the references from 
which these were drawn; the table includes several sources 
that were not investigated in IWGP13I because they have no 
PRE bursts in the MINBAR catalog. 

2.1. Burst selection 

We restrict our sample of type I bursts to exclude events that 
are unsuitable for analysis. We discarded 57 bursts for which 
no Standard-2 data was available, preventing estimation of the 
instrumental background. Some burst sources lie in crowded 
fields containing other LMXBs within the ~ 1 0 RXTE field of 
view. If the other source(s) were active at the time of obser¬ 
vation, then their persistent emission would be confused with 
that of the burst source and it would be impossible to separate 
the persistent emission of the burst source from that of the 
other source(s). Such bursts need to be excluded from consid¬ 
eration. Our procedure for removing source-confused bursts 
is the same as IWGP13I §4.1. A total of 168 bursts were ex¬ 
cluded as being source confused. Our catalog also contains 16 
type I bursts from the Rapid Burster taken during offset point¬ 
ings to avoid confusion with the nearby 4U 1728-34. We also 
excluded 35 bursts whose radius expansion status was unclas- 
sifiable due to data gaps or very low flux, and 58 bursts clas¬ 
sified as "marginal" (i.e. satisfying only some of the criteria 
of §2.3 inlCXM). 

2.2. Characterizing the persistent emission 

As in the previous studies lG08l and IWGP13I we adopted 
the integrated X-ray flux for a 16-second interval prior to 
the start of each burst as the persistent emission. This spec¬ 
trum includes a time-dependent contribution estimated as in 
IWGP13L Subsequent model fits to each persistent (and burst) 
spectrum used the corresponding model spectrum estimated 
for that burst as background. 

For each burst we fit the persistent emission with a set 
of nine alternative models in turn. These are summarized 
in Table [2] and are the same six persistent models as in 
IWGP13L with the addition of wabs*disko, wabs*diskm, 
and wabs* (CompTT+gaussian). The first two of these 
include accretion rate explicitly as a variable parameter, mak¬ 
ing it possible to investigate the relationship between the nor¬ 
malization factor f a and M. The latter was used bv lPeille et al.l 
Uni and we include it for completeness. We then selected 
the fit that gave the best (i.e. lowest) xt to represent the per¬ 
sistent spectrum for that burst. As shown in Figure [Q this 
suite of models provide statistically adequate fits so we re- 

4 see http://burst.sci.monash.edu/minbar 
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TABLE 1 

N h Values and Eddington Fluxes 


Source 

Non-PRE Bursts 

PRE bursts 

SWFC a 

SBC b 

Nh 

10 22 cm^ 2 

Reference 

^Edd 

10 -9 erg s -1 cm -2 

4U 1746-37 

22 

3 

0 

2 

0.3 

i 

5.8 ± 1.4 

MXB 1659-298 

11 

12 

0 

0 

0.2 

2 

17.0 ±4.0 

Aql X-l 

52 

12 

0 

0 

0.4 

3 

99.6 ±21.3 

SAX J1810.8-2609 

4 

1 

1 

0 

0.3 

4 

111.2 ±2.5 

IGRJ17480-2446 

277 

0 

0 

0 

0.5 

1 

C 

CygX-2 

34 

7 

0 

0 

0.1 

5 

13.1 ±2.1 

GS 1826-24 

67 

0 

0 

0 

0.4 

6 


XTEJ1709-267 

2 

0 

0 

0 

0.4 

7 


4U 1820-30 

0 

16 

34 

7 

0.2 

1 

57.1 ±7.7 

SAX 1808.4-3658 

1 

8 

2 

0 

0.1 

8 

230.1 ± 13.2 

GX 17+2 

5 

2 

0 

4 

1.9 

9 

13.3 ±2.5 

Rapid Burster 

19 

0 

0 

0 

1.7 

10 


IGR J17473-2721 

2 

0 

0 

0 

3.8 

11 

113.5 ± 12.1 

Cir X-l 

13 

0 

0 

0 

0.7 

12 


HETEJ1900.1-2455 

0 

7 

0 

0 

0.2 

13 

123.9 ± 10.6 

4U 2129+12 

5 

1 

1 

1 

0.0 

1 

40.8 ± 1.6 

1M 0836-425 

17 

0 

0 

0 

2.2 

14 


4U 1735-444 

7 

10 

0 

0 

0.1 

15 

34.2 ± 5.6 

1A 1744-361 

3 

0 

0 

0 

4.5 

16 


4U 1722-30 

0 

3 

23 

1 

0.8 

1 

61.7 ± 12.4 

GX 3+1 

1 

2 

0 

1 

1.6 

17 

59.9 ± 0.9 

EXO 1745-248 

15 

2 

0 

0 

3.8 

1 

69.0 ± 0.2 

XTEJ 1739-285 

6 

0 

0 

0 

2.0 

18 


2S 0918-549 

1 

2 

3 

0 

0.3 

19 

119.2 ± 12.4 

Ser X-l 

11 

6 

0 

0 

0.4 

5 

29.4 ±7.1 

GRS 1741.9-2853 

1 

4 

1 

0 

11.3 

20 

35.3 ± 10.9 

2E 1742.9-2929 

65 

1 

0 

0 

1.2 

18 

37.8 ± 1.4 

4U 1728-34 

30 

64 

0 

0 

2.6 

21 

95.0 ± 8.4 

XB 1832-330 

0 

1 

0 

0 

0.1 

1 

33.7 ±4.4 

IGRJ 17511-3057 

7 

0 

0 

0 

0.6 

22 


4U 1608-522 

28 

18 

5 

0 

0.9 

23 

167.2 ± 26.0 

GRS 1747-312 

1 

1 

0 

0 

1.4 

1 

13.4 ±4.4 

XTE 1814-338 

27 

0 

0 

0 

0.2 

18 


XTE J1723-376 

2 

0 

0 

0 

7.9 

18 


4U 1254-69 

3 

0 

0 

0 

0.3 

24 


4U 1702-429 

43 

5 

3 

0 

1.9 

18 

87.7 ±4.5 

SAX J1750.8-2900 

3 

2 

1 

0 

0.9 

18 

54.1 ±2.1 

XTEJ 1701-462 

1 

2 

0 

0 

2.0 

25 

43.4 ± 1.4 

XTE J1810-189 

3 

1 

0 

0 

4.2 

26 

54.2 ± 1.8 

EXO 0748-676 

143 

5 

0 

0 

0.8 

27 

46.5 ±4.6 

4U 1916-053 

1 

12 

0 

0 

0.3 

28 

30.6 ± 3.6 

XTE 2123-058 

3 

0 

0 

0 

0.1 

29 


XTE 1759-220 

5 

3 

0 

0 

2.8 

18 

15.7 ±0.8 

4U 1636-536 

269 

76 

2 

0 

0.2 

5 

72.6 ±9.1 

IGR 17191-2821 

5 

0 

0 

0 

0.3 

30 


SLX 1744-300 

10 

0 

3 

0 

4.5 

31 

13.9 ±3.1 

SAX J1748.9-2021 

16 

11 

0 

0 

0.5 

1 

38.0 ± 6.0 

4U 1323-62 

35 

0 

0 

0 

2.4 

32 


KS 1731-260 

21 

4 

3 

0 

1.3 

33 

48.6 ±5.6 

XTE J1710-281 

37 

3 

0 

0 

0.4 

34 

7.1 ± 1.5 

IGRJ 17498-2921 

0 

1 

0 

0 

1.2 

18 

51.6 ± 1.6 

SAX J1747.0-2853 

6 

10 

2 

0 

8.8 

35 

52.5 ±7.1 

SLX 1735-269 

1 

0 

0 

3 

0.1 

36 

49.4 ± 3.9 

4U 1705-44 

77 

4 

0 

0 

1.9 

37 

41.0 ±3.8 

4U 0513-40 

11 

4 

3 

0 

0.0 

1 

14.5 ±3.5 

SAX J 1806.5-2215 

4 

0 

0 

0 

1.0 

18 
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Fig. 1.— Distribution of for the spectral fits to the persistent emission 
(histogram) compared to a theoretical distribution of xt f° r a collection of 
spectra with the same number of degrees of freedom (curve). A Kolmogorov- 
Smirnov test gave a 97% probability (D = 0.11 (that the two curves are con¬ 
sistent with having been drawn from the same distribution, indicating that 
our suite of models for the persistent emission spectra is adequate for use in 
subsequent work. 

gard these persistent models acceptable in the analysis of the 
bursts themselves. 

There is no difference in the distributions of best-fitting per¬ 
sistent models between the PRE and non-PRE bursts. We as¬ 
signed each model a distinct numerical label and performed a 
Kolmogorov-Smirnov test on the resulting distributions, with 
D = 0.03 and a 95% probability that both are consistent with 
having been drawn from the same distribution. 

A gas pressure dominated accretion disc (wabs *diskm) is 
not the preferred model for any persistent emission spectrum, 
indicating that this is not a good description of accretion discs 
in LMXB systems. Instead, we expect the disc to be radiation 
dominated, and possibly showing hard emission. 

For 6 bursts, no persistent emission model could be fit with 
xl < 3.5. These bursts were excluded from further consider¬ 
ation. 

2.3. Fitting procedure 

As in lWGP13l we first fit every burst spectrum with the stan¬ 
dard approach: a blackbody spectrum, corrected for interstel¬ 
lar absorption using the appropriate hydrogen column density 
for the burst source (see Table [T]). We use the recorded pre¬ 
burst emission, including the instrumental background, as the 
background for these initial fits. See IG08I for further details 
regarding this fitting method. The standard fits provide a com¬ 
parison model and spectral parameters, as well as sensible 
spectral parameters to seed the variable persistent emission 
fits. 

We then re-fit all the burst spectra, replacing the back¬ 
ground with the instrumental background, and adding a vari- 
a ble p ersistent emission component appropriate for that burst 
( 32,21) . That is, we modelled the total spectrum as 

S(E) = A(E) x B(E-,T bh ,K bh ) + f a x P(E) + b(E) inst , (1) 

where S(E) is the fitted spectrum as a function of energy E, A 
is the absorption correction due to interstellar hydrogen (see 
uh values listed in Table [T]). B is a blackbody spectrum with 


temperature T bb and normalization K bb , and b ms i is the instru¬ 
mental background. P is a model for the pre-burst persistent 
emission that also includes absorption, though for some per¬ 
sistent models we do not retain the same nn~ this is a means 
to adjust the low energy end of the accretion spectrum and 
should not be interpreted as saying anything physical about 
the hydrogen column density. The accretion enhancement is 
quantified with the f a normalization factor. We have fit a total 
of 160,017 spectra. 

3. RESULTS 

Profiles of fit parameters, using both methods, for three ex¬ 
ample bursts are shown in Figure [2] We have selected bursts 
from three different sources to illustrate the results. These 
plots clearly show that f a is enhanced to several times the pre¬ 
burst level, as found for radius expansion bursts by 1WGP13L 
This result is typical for all the bursts in our analysis. The 
f a for spectra preceding the burst appear to be slightly lower 
than unity, but this is an artefact of fitting these spectra with a 
model that includes a burst component which in reality is ab¬ 
sent before the onset of nuclear burning. This causes some of 
the persistent flux to be misidentified as burst emission. Fit¬ 
ting the pre-burst spectra with just a normalization-variable 
persistent model gives results consistent with unity. 

In Figure [3] we show the distributions of xt for the variable 
persistent normalization fits and the standard approach fits, 
compared with a theoretical distribution of xt for a model 
that adequately describes the data. Fits from both PRE bursts 
and non-PRE bursts are shown. We performed Kolmogorov- 
Smirnov tests on the measured distributions of xt against the 
theoretical. D and p values are listed on Figure 0 with the 
variable fits listed first. These results confirm what is visually 
evident in the Figure: the variable persistent normalization 
fits significantly improve the quality of the spectral fits for 
both radius expansion and non radius-expansion bursts. The 
spectral fits for PRE bursts are generally poorer than non-PRE 
bursts for both fitting methods and, although the variable per¬ 
sistent normalization method improves the distribution of xt> 
the deviation from a model that is statistically consistent with 
the data is still present. It is clear that radius expansion intro¬ 
duces a significant spectral effect on top of the variations in¬ 
duced by enhanced persistent emission. The second and third 
panels show that cooling tail spectra of PRE bursts are more 
poorly fit than spectra from non-PRE bursts. This may imply 
that PRE affects the disc and/or photosphere in a manner that 
persists for some time after the photosphere has returned to 
the surface of the star. 

We assessed the relative improvement of the variable per¬ 
sistent flux treatment over the standard analysis by estimating 
the Bayes factor, as follows. We interpret the probabilities 
calculated from the K-S statistics as the global likelihood for 
each model, i.e. p(D\M,,I) where D is the set of x_ 2 values 
corresponding to the model fits, M, represents each fitting ap¬ 
proach (model), and I the priors. Even if this interpretation 
is not strictly correct, we expect that the K-S probability is 
proportional to the actual likelihood, so that the ratio of prob- 
abi lities is a good estimate of the ratio of likelihoods. Follow¬ 
ing [Gregory 0 200 51), the other contribution to the Bayes factor 
calculation is the penalty imposed by the additional degree 
of freedom introduced with the variable persistent flux factor 
f a . Since the standard approach is a special case of the vari¬ 
able persistent flux approach (with f a fixed implicitly at 1), 
we compare the effective uncertainty for the former approach 
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TABLE 2 

XSPEC MODELS FOR FITTING THE PERSISTENT EMISSION 


XSPEC model 

Number 

Notes 


of spectra 


wabs*disko a 

160 


wabs*diskm a 

0 


wabs*(bbodyrad+powerlaw) 

425 

nH fixed to literature values 

wabs*(bbodyrad+powerlaw+gaussian) 

271 

nH fixed, Gaussian energy set to 6.4keV 

wabs* (compTT b ) 

39 

nH allowed to vary 

wabs*(bbodyrad+powerlaw) 

46 

nH allowed to vary 

wabs*(gauss+bremss) 

321 

all parameters variable 0 

wabs*(bbodyrad+diskbb d ) 

183 

all parameters variable 0 

wabs*(compTT+gaussian) e 

314 

nH fixed, Gaussian energy set to 6.4keV 

Total usable spectra 

1759 


Rejected due to source confusion 

res 

Other active sources in field 

Background data missing or unusable 

57 


No good persistent model fit 

6 

Minimum xj >3.5 

Marginal radius expansion 

58 


Radius expansion status unclassified 

35 



a SeeiSt ella & Rosneii il 984 i) 
b See lTitarchukjfl99l 

c All parameters are variable for the generation of persistent emission models. Their values are subsequently frozen for the 
burst spectral fits in ClU 

d See XSPEC manual ihttp: //heasarc. gsfc . nasa. gov/xanadu/xspec/manual/manual. html'i and refer¬ 
ences therein 
c See lPeille et aTH20l3) 


(arising from the Gaussian statistics on the pre-burst persis¬ 
tent spectra) with the prior for the latter approach, which is 
flat between f a = -100,100. The median total counts in all 
the pre-burst spectra is approximately 9,000 counts, and so 
we adopt 10 4 counts, giving a characteristic width of the prior 
in the standard approach of 5f a /f a = v 7 10 4 /10 4 = 0.01. The 
overall estimate of the Bayes factor is then 

p(D\M u I) 5f a _ 5.77 xlO- 3 0.01 
~ p(D\M 1 J) A/ a 4.51 xl0“ 9 200 1 

Thus, even taking into account the additional degree of free¬ 
dom introduced by allowing the normalisation of the persis¬ 
tent flux to vary, we conclude that the variable persistent flux 
improves on the standard approach with a Bayes factor of 64. 

The addition of a third variable parameter to the model can 
have an adverse effect on the uncertainties in the fit parame¬ 
ters. In Figure [I] we compare the uncertainties in fit parame¬ 
ters between the two approaches. It is clear that introducing f a 
has a slight adverse effect on the uncertainty in the blackbody 
temperature. The uncertainties in blackbody temperature are, 
on average, 25% larger for the variable persistent flux method. 
The uncertainty in the blackbody normalization on average in¬ 
creases by 49%, to about 18% of the measured value. The un¬ 
certainty in f a is on average 27% of the measured value. The 
poor constraints on f a and blackbody normalization are likely 
due to th e spe ctral similarity of the two components. Recently 
iBarriere et al.l (120141) showed a severe example of this prob¬ 
lem in using the variable persistent flux method in a NuSTAR 
observation of GRS 1741.9-2853, in which the f a factor was 
not constrained at all. As they point out, the problem occurs 
when the pre-burst persistent emission is faint, and when it 
closely resembles the burst component in spectral shape. 

IWGP131 found, for PRE bursts, only a weak correlation be¬ 
tween f a and burst flux. They found that, although most of the 
high f a values occurred during radius expansion, there was 
nonetheless great variability in f a when the burst flux was rel¬ 
atively constant. We have repeated this analysis for non-PRE 
bursts. We have taken all spectra for each burst between the 


beginning of the burst and the point at which burst flux de¬ 
clines below 10% of its peak value, and performed Kendall 
r rank correlations on burst flux against f a for each burst. 
We found positive correlations in 1000 out of 1419 non-PRE 
bursts, of which 919 are at greater than 3er significance. Per¬ 
forming the same test on the cooling tails of PRE bursts, we 
found positive correlations in 299 out of 316. Of these, 286 
were at greater than 3cr significance. These numbers demon¬ 
strate that there is a more direct relationship between burst 
fl ux a nd f u for non-PRE bursts than for PRE bursts. 

IWGP13I found that, in PRE bursts, the variations in mea¬ 
sured f a cannot be attributed to counting statistics. We re¬ 
peated this test for all non-PRE bursts in our catalog. For 
all spectra obtained during non-PRE bursts we generated 100 
simulated spectra using the standard approach spectral param¬ 
eters, and incorporating counting statistics typical of the de¬ 
tector. We then fit these simulated spectra with the variable 
persistent normalization method, obtaining a standard devia¬ 
tion of f a measurements arising from counting statistics in the 
absence of any variation of persistent flux. We found real f a 
measurements greater than unity to 5er significance in 214 of 
62783. In the absence of a real effect we would expect at most 
one such case arising by chance. This experiment shows that 
enhanced persistent emission is present in non-PRE bursts as 
well. 

Since PRE bursts are by nature brighter than non-PRE 
bursts, it is possible that the difference in xt between the 
two classes is simply due to superior signal-to-noise in the 
brighter bursts. To investigate this we took a subset of spectra 
from both PRE and non-PRE bursts, with bolometric fluxes 
between 4.0 x 10~ 8 and 5.0 x 10 -8 erg s -1 cm -2 and blackbody 
temperatures between 2.0 and 2.5 keV, these quantities being 
measured with the standard approach. There were 331 spectra 
from PRE bursts and 655 spectra from non-PRE bursts in this 
segment of the parameter space. These had mean % 2 of 1.17 
and 1.35 for non-PRE and PRE bursts respectively, and a K-S 
test indicated a probability of 1.79 x 10 -6 that these % 2 are 
drawn from the same distribution. This shows that the poorer 
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Fig. 2.— Comparison of fitting a variable persistent emission factor, f a , to the standard approach fits for non-PRE bursts from three different neutron stars. 
We have selected a typical hydrogen accretor (4U 1608-52), a dipping source (EXO 0748-676), and a He-accretor and likely ultracompact binary (4U 1728—34). 
The variable fit approach yields consistently lower unabsorbed burst component fluxes (top panels). The contribution to the flux from the variable persistent 
emission increases to several times the pre-burst level (middle panels) during the rise of the burst. Allowing the persistent emission to vary improves the xt 
(lower panels). 


spectral fits generally observed in PRE bursts are not due to 
better signal-to-noise in the bright bursts. 

4. IS THE PERSISTENT EMISSION SPECTRAL SHAPE CONSTANT? 

Our approach of treating the persistent emission as a 
contribution of fixed shape depends upon the assumption 
that its shape does not change, at least not on timescales 
of the order of the burst duration, _or to a degree suffi¬ 
cient to affect our spectral fits. IWGP131 did not investi¬ 
gate this matter directly, but their finding that f a generally 
returns to its pre-burst level, with good sugg ests that 
this as sum ption is not unreas onable. Similarly, iKeek et all 
(12014ft and iPeille et all (12014ft also report that the persistent 
emission returns to its ori ginal level soon afte r the burst. 
iThomnson et all (12005ft and fBagnoli et all (12013ft have found 
that the persistent emissions of 4U 1826-24 and the Rapid 
Burster respectively are remarkabl y stab le over kilosecond 
timescales between bursts. iLinares et al.l (12014ft found vari¬ 
ability over ~10 minute periods in the intensity of the per¬ 
sistent emission of IGR J18245-2452, with no change in the 


shape of the persistent emission. Dipping systems, such as 
EXO 0748-676 and X 1658-298, show somewhat more vari¬ 
ability (ISidoli et al Jl2005t lOosterbroek et al.1l2001bl) on~l m 
timescales, and the spectrum becomes harder during dips, as 
indicated by hardness ratios. However, the dips are due to 
obscuration in these high-inclination systems. Outside of dip¬ 
ping intervals, the hardness ratios are effectively constant. It is 
therefore likely that no change in the shape or intensity of the 
persistent emission is to be expected in the absence of burst 
luminosity on minute to hour timescales. 

On shorter timescales, the pre-burst light curvejs remark¬ 
ably flat preceding b ursts fo r 4 U 16 36-536 (iN'ath et al.l2002l) . 
the Rapid B urster dBagnoli et al.l 120131) . and 4U 1608-52 
(IPeille et al.ll2014ft . and is also flat over a large rang e of ac¬ 
cretion rates in IGR J17480-2446 (iMotta et al. 11201 ll) . These 
works suggest that short-term intensity variations are not 
prevalent in the persistent emission preceding bursts. To gain 
further confidence that the persistent emission shape and in¬ 
tensity are not significantly variable on shorter timescales than 
minutes or hours, in the absence of a burst, we took every 
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Fig. 3.— Histograms of the xl> f° r the variable persistent normalization and standard approach spectral fits (red and black histograms respectively), and the 
theoretical distributions for a model that adequately describes the data (red and black curves respectively). The first panel shows spectra from Eddington- 
limited PRE bursts, the second panel shows cooling tail spectra from PRE bursts, and the third panel shows spectra from non-PRE bursts. It is obvious that the 
variable persistent normalization approach gives better fits overall than the standard approach. It can also be seen that a blackbody burst spectrum is far from 
acceptable in PRE bursts, even in the cooling tail when the atmosphere has returned to the surface of the star. The variable persistent normalization approach 
corrects for much of this effect. Non-PRE spectra are better described by a blackbody and constant persistent model, but even these are significantly improved 
with the variable persistent normalization approach. Results of Kolmogorov-Smirnov tests are also shown; we list the K-S statistic D and null hypothesis 
probability p for the variable persistent normalization and standard fits respectively. 




Fig. 4.— Histograms of the relative uncertainties in the spectral parameters using both approaches. The addition of the f a parameter causes the temperature 
and normalization of the blackbody component to be less tightly constrained than in the standard approach. The persistent flux factor f a itself is often very poorly 
constrained. 
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Fig. 5.— Measured distribution of /„ (red histogram), and measured and 
theoretical Xu (thick and thin blue curves) for pre-burst emission in the ab¬ 
sence of any nuclear burning. The /„ values are very strongly peaked around 
1. and the measured xt a| c low. indicating that the pre-burst emission does 
not change intensity or spectral shape to any detectable degree. 

0.25s pre-burst spectrum from every burst and fit it with just 
f a times the persistent model for that burst, i.e., no thermonu¬ 
clear burst emission. The distributions of f a and xf are shown 
in Figure 0 This data typically covers a few tens of seconds 
prior to every burst. It is plain to see that f a is distributed 
tightly around unity, indicating that the persistent emission 
generally does not vary greatly in intensity on timescales of 
1 s or less. The xt distribution also implies that the shape of 
the persistent emission is not changing appreciably, although 
of course the short exposure times and associated low photon 
counts make detecting such changes difficult. 

To test for more subtle spectral shape changes, we increased 
the exposure time to obtain better signal-to-noise. For the 
three bur sts stu died i n Figu re E and the fifteen bursts investi¬ 
gated bv lPeille et al.l (120141) (see their Table 1), we inspected 
the RXTE Standard-2 data covering the entire observation. 
These are spectra binned into 16s intervals. We fit each of 
these spectra with the persistent model for that burst, with all 
parameters frozen, and f a normalization allowed to vary. We 
considered only spectra recorded within an hour before or af¬ 
ter the burst but excluding those during the burst, to avoid 
contamination by the burst emission. That is, we excluded 
spectra recorded between the beginning of the burst (burst 
flux has reached 1/4 of the peak burst flux) and the end of 
the PCA data of the lG08l catalog. This data typically covers a 




Fig. 6.— Measured distribution of /„ (red histogram), and measured and 
theoretical xt (thick and thin blue curves) for Standard-2 data covering the 
observations of eighteen selected bursts. The f a values are tightly clustered 
around unity, and the xt are marginally higher than the theoretical distribu¬ 
tion assuming a well-fitting model. This suggests that the persistent emission 
does not change shape or intensity to any detectable degree on timescales of 
an hour before and after a type I burst. 

few minutes to an hour before and after the burst. We also ex¬ 
cluded an eclipse of EXO 0748-676 beginning at 2004-09-25 
07:20:29, 15 minutes before the burst, when no flux was ob¬ 
served from the neutron star (see also lHoman et al.l2003l) . and 
spectra occurring after a data gap towards the end of an ob¬ 
servation of 4U 1636-536 that began at 2001-09-05 04:56:51. 
We analyzed the distribution of f a and xt- These are plotted 
in Figure [6] showing that almost all of f a values measured 
this way are clustered tightly around f a = 0.95 - 1.0. There is 
a smaller peak at around f a ss 0.85, dominated by two bursts, 
but visual inspection of their f a curves show their f a values to 
stay approximately constant throughout the observation. 

The steadiness of both f a and xl over time suggests that the 
persistent emission does not change greatly in either shape or 
intensity on time scales up to an hour before and after a type 
I burst. 

Finally, we investigated the extent to which we can de¬ 
tect changes in the accretion spectral shape if they really are 
present, assuming that a change in the intensity of the per¬ 
sistent emission reflects a change in the accretion rate. This 
is particularly necessary near the burst peak, since this is 
where the persistent emission spectral shape is most uncer¬ 
tain; we have already found that it returns to its original shape 
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deep in the cooling tail. We selected a non- PRE burst from 
4U 1636-536 that uses the disko model (IStella & Rosnerl 
fl984l) . recorded 2005 Apr 4, 10:48:43. For this burst PCUs 
0 and 2 on RXTE were active. This model contains M/Meaa 
explicitly as a variable parameter and incorporates physically- 
motivated changes in the spectral shape. This model is as¬ 
sumed to be in equilibrium for this burst and does not take 
into account changes in its structure induced by a burst; there 
are as yet no theoretical spectra of a disk modified in this way, 
and we expect that the associated spectral shape will deviate 
from a blackbody to a similar degree, though not perhaps in 
the same manner. 

The best fit accretion rate was 0.205M/MEdd and the peak 
f a was 3.3. The normalization of this persistent model was 
2 cos i/d 1 = 1.648, where i is the inclination and d is the dis¬ 
tance in units_of 10 kpc. At an assumed distance of 6 kpc 
dGallowav et all 120061). this gives an i nclination of 72°, i n 
agreement with fPandel et all (120081) and lCasares et all (120061) . 
who also infer a high inclination (greater than 64°, and 
36° -74° respectively) for this system. The normalization was 
kept frozen in the subsequent analysis. We generated 2,000 
simulated spectra, folded through a PCA response appropri¬ 
ate to the observed burst, of the wabs* disko model withM 
taking random values between 0.067MEdd and 0.667MEdd- We 
then fitted those simulated spectra with f a times the original 
persistent model. The measured f a and xt are given in Fig- 
ure[7] left two panels. In the absence of burst emission we can 
clearly detect a change in spectral shape via the increase in xj, 
with increasing M as well as an increase in the intensity of the 
persistent emission. We then repeated the process, this time 
including a blackbody with temperature 2.0 keV and normal¬ 
ization representing a sphere with radius 10 km at a distance 
of 6 kpc (i.e., the surface of 4U 1636-536). This we fit with 
the usual variable persistent flux approach (right two panels). 
It is evident that, although we can still detect enhanced persis¬ 
tent emission through an increased f a , we have lost the ability 
to detect a change in its spectral shape. This experiment gives 
us confidence that our approach of simply varying the normal¬ 
ization of the persistent emission will introduce no detectable 
systematic effects. 

This analysis is not sensitive to changes in the persistent 
spectrum o utsi de the 2.5-2 0 keV energy range we consider in 
this paper, lin’t Zand et all (2013) combined data from Chan¬ 
dra, which is sensitive at energies below 2keV, with RXTE 
data for a radius expansion burst from_SAX J1808.43658. 
They obtained similar f a values as IWGP13I did studying the 
same burst with only RXTE data, suggesting that changes in 
the persistent spectrum below 2.5keV are negligible. At high 
energies, there is evidence that the hard X-ray (>30 keV) flux 
decreases during a burst, in apparent c onflict with our result 
(IChen et aT]l2013t [ji et al.ll20T3Ll2014aIB) . However, this phe¬ 
nomenon is attributed to the rapid cooling of the accretion 
disc corona, which contributes the majority of the very hard 
X-rays. An increase in 2.5-20keV flux attributed to accretion 
rate change does not conflict with a simultaneous quenching 
of high energy photons through the cooling of the corona. 

5. EDDINGTON FLUXES FROM NON-PRE BURSTS 

For neutron stars for which no radius expansion bursts 
have been observed, the only lower bound on their Edding¬ 
ton fluxes are the peak fluxes of their brightest bursts. A 
better knowledge of their Eddington fluxes would be of ob¬ 
servational utility, as many of these sources show interesting 


properties: the prolific burster GS 1826-24 exhibits very reg¬ 
ular bursts that are in close agreement w i th numerical models 
of X- ray bursts (e.g.. iHeger et al.l 120071 : iGallowav & Lampel 
HHJ). The 11 Hz pulsar IGR 17480-2446 is t he most slowly 
r otati ng source that shows burst oscillations (ICavecchi et alJ 
1201 ll) . making that source an important test of burst oscilla¬ 
tion models. Neither source has a single reported PRE burst. 
The Rapid Burster (MXB 1730-335) is our major source of 
knowledg e about type II bursts, and has one reported PRE 
burst (ISala et al.ll2012h . This event, however, was detected by 
the Swift X-ray instrument and the measurement was not sen¬ 
sitive enough to constrain the Eddington flux to better than 
about 50% (see their Fig. 2). 

An interesting observation is that all the results of lWGP13l 
fall below the line f a 7 = 1 , where 7 is the pre-burst accretion 
flux as a f raction of the Eddingt on fl ux, consistent w i th pre - 
dictions of lBurger & Katzl(Il983h and lMiller & Lambl(Il996h : 
both works suggest that the Eddington accretion rate is a nat¬ 
ural upper limit to M. This observation raises the possibility 
of obtaining a lower bound on the Eddington flux for a given 
source, that is, /Aid > f a F peK , where Tp ers is the bolometric 
flux of the pre-burst persistent emission. 

In Figure[ 8 ]we show peak f a against 7 for every type I burst 
detected by RXTE. For 0.01 < 7 < 0.4 the non-PRE bursts ap¬ 
pear to be bounded above by f a 7 < (. f a i)max , where (J a l)max 
appears to be somewhat less than 0.5. This would give a lower 
limit on the Eddington flux roughly twice as constraining as 
the most conservative estimate f a F pers < FecIcI- Determining 
7 requires knowing /Aid; our method for determining the Ed¬ 
dington fluxes for PRE bursters is given in the Appendix. 

We now describe a method by which we estimate the upper 
bound of f a 7 that is reproducible and more objective than es¬ 
timating by eye. We treat the measured x; = (f a ^)i values as 
normal distributions centered around their measured values, 
with standard deviations 07 equal to the uncertainty in the /',, 
times 7. This normal distribution is given by 


G{x,(j,.Xi) = 


1 

~ff = ex P 

y lira i 


( x-Xj ) 2 

2(7; 


(3) 


The fraction of this distribution lying above x max = {f,rf)n 


A((7;,X;) = / G(x, (7;,X;)(/x 

x ; X/h/ix 


1 +erf 


1/2(7; 

and so the total over all data points is 


1 


AT = ,£l + erf 


Xf-x, 


1 •*'max 


. \/2<Ti J ’ 


(4) 


(5) 


where n is the number of data points. 

The question is how to find x max . First we make a guess ( as 
to the value of x max and calculate an At from this and the data 
points X,. Then we take a set of n data points x, with random 
values between 0 and £. We assign to each x; an uncertainty 
of (7,-, that is, a random point is given a measured uncertainty. 
This assumes the real data points are evenly distributed with 
uncertainties independent of their values, but the procedure 
can be modified to suit other assumptions. From this random 
set x; we calculate an expected total area (A j(rand)) that can 
be compared to the real one (At (observed)). If, for instance. 
At (rand) < At {observed) then our guess for C, was too high 
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Fig. 7.— Results of fits to simulated disko spectra with M allowed to vary, both with (panels a. and b.) and without (panels c. and d.) a burst component. In 
both cases an increase in f a is detected. When there is no burst component, and only a change in the intensity of the persistent emission is being tested for, the 
increasing xi with increasing M indicates that changes in the spectral shape are detectable. In the presence of burst emission, it is no longer possible to discern 
the spectral shape changes in the persistent emission. We can therefore simply vary the normalization of the persistent emission without introducing systematic 
effects. 


and we try again with a lower £. Once the £s converge, we 
take this value to be x max . Making many realizations of the 
random set, we can get an expected x max . However, to avoid 
the analysis assigning too much weight to outliers, we take 
many bootstrap subsamples of the original data set and take 
the overall x max to be the average of those from the bootstrap 
samples; it also gives us uncertainties on x max . 

Applying this procedure, we get x max = 0.84 ± 0.34 for the 
PRE bursts and x max = 0.37 ± 0.089 for the non-PRE bursts. 
We note that the f a 7 relation is a phenomenological one and 
does not depend on any physical interpretation, though it is 
consistent with the requirement that the accretion luminos¬ 
ity not exceed Eddington. To test whether our new method 
gives a more stringent lower bound than the peak burst flux, 
we compared the two methods for every source for which ten 
or more non-PRE bursts have been observed, but excluding 
the rapid accretors Cyg X—2 and GX+17-2. We also exclude 
SLX 1744-300, as that source is difficult to distinguish from 
the nearby (A 9 « 3') SLX 1744-299 (ISkinner et al.lfl990h . 
We first apply the method to sources for which the Eddington 
flux is known, to make sure the new method never overesti¬ 
mates the Eddington flux. In 2 of 16 sources, the new method 
gave a greater lower bound than the brightest non-PRE burst, 
and in no case did the new method overestimate the Edding¬ 
ton flux beyond observational uncertainties. These results are 
shown in Figure [9] 

Applying the new method (with x, = (/ fl Tp ers ),- since 7 
is unknown) to sources with unknown Eddington flux, we 
get a better bound than the non-PRE bursts for the sources 
Cir X—1 and IGR 17480-2446, with /Aid greater than 


18.1 ±5.5 x 10 “ 9 and 42.0± 10.5 x 10 “ 9 erg s ” 1 cm " 2 re¬ 
spectively. IGR 17480-2446 is located in the globular 
cluster Terzan 5, with a known distance of about 5.5 kpc 
(e.g. JPapitto et al.l 12011 ). and a hydrogen accreting neutron 
star with 1.4M 0 would have an Eddington flux of 57 x 10 -9 
erg s " 1 cm~ 2 . Cir X-l is b elieved to be 7.8-10.5 kpc distant 
dJonker & Nelemansil2004l) . Our result suggests it is probably 
at the nearer end of that range, since the more distant value 
would give an Eddington flux of 16 x 10 ~ 9 erg s _I cm -2 for a 
canonical neutron star, slightly higher than our estimate. We 
also get a marginally more constraining Eddington flux for 
1M 0836-425 of 23.1 ±6.5 x lO^ergs " 1 cm" 2 . 

6. CONCLUSION 

We h ave extended the vari able persistent flux method devel- 
oped bv lWorpel et al.l (120131) and lin’t Zand et all (120131) to all 
type I bursts observed by RXTE. Our method gives superior 
spectral fits to type I burst spectra, whether these are radius 
expansion bursts or not. Our detailed conclusions are as fol¬ 
lows: 

i) The variab le persistent normalization approach devel¬ 
oped in lWGP13l applies as well to non-PRE bursts as it does 
to radius expansion bursts. The quality of spectral fits gen¬ 
erally improves, and the intensity of the persistent emission 
usually increases during a burst, typically by a factor of 3-4. 

ii) The addition of a third variable parameter increases the 
uncertainties in the fit parameters. In particular there is a 
degeneracy between f a and blackbody burst component nor¬ 
malization, likely due to the indistinctness of the two spectral 
components. 
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Fig. 8 .— Peak f a against 7 for every burst in the catalogue from sources for which Eddington fluxes are listed in |WGPl"3l PRE bursts are shown in red, and 
non-PRE bursts are shown in blue. There is a clear division between the two classes of bursts. PRE bursts are bounded above by f a 7 < 0.84, while the non-PRE 
bursts are bounded above by fa'y < 0.37 (as calculated by the method explain ed in ©■ The dotted lines indicate 1 rj intervals estimated using the bootstrap 
method. We also show the hard upper bound fa'y < 1 (e.g. IBurger & Katzi 19831) . The clump of points at 7 Ri 1 represents bursts from the rapid accretors Cyg X-2 
and GX 17+2; these are atypical bursters and we do not consider their bursts in deriving fa'y relations. 


iii) Photospheric radius expansion adds a confounding 
spectral effect, likely caused by the changing structure of the 
photosphere. When this is corrected for, by considering only 
non-PRE bursts, the spectral fits are much improved, to a 
Bayes factor of 64. 

iv) Any remaining discrepancy from a blackbody plus vari¬ 
able persistent normalization model is likely to be due to a 
slight non-blackbody character to the burst emission, rather 
than spectral shape changes in the persistent emission. This is 
because changes in the shape of the persistent emission are not 
detectable in the presence of an additional burst component 
with the quality of the data currently available. Investigations 
into the response of the accretion spectrum to a burst tran¬ 
sient are desirable, but must await more sensitive ins trum ents 
suc h as the Large Observatory for X-ray Timing dFeroci et al.l 

I20l2h . 

v) We have found no evidence that the persistent emission 
spectrum changes shape on timescales of up to an hour before 
or after a burst, and no major changes in its normalization. 

vi) For sources that have never shown a radius expansion 
burst, we can obtain a lower bound on their Eddington lumi¬ 
nosity from the observation that the persistent emission en¬ 
hancement for non-PRE bursts seldom exceeds a certain frac¬ 
tion (about 37%) of the Eddington flux. For three sources, 
this new lower bound was found to be a better constraint than 
the peak fluxes of the brightest (non-PRE) bursts for those 
sources. 

We recommend that varying the normalization of the pre- 


sistent flux be adopted as standard for analysing type I bursts, 
though care must be taken to account for the increased uncer¬ 
tainties in the fit parameters caused by the introduction of a 
third parameter. 

APPENDIX 

Table Q] lists Eddington fluxes for all sources from which 
at least one PRE burst has been observed by RXTE’s PCA 
or WFC instruments. This data is supplemented by a small 
number of PRE bursts recorded by other X-ray satellites. 
These are mostly medium to long duration bursts and are 
listed in Table 0 This sample is heterogeneous, but repre¬ 
sents our only knowledge of the Eddington flux of the source 
SLX 1735-269, where only non-PRE bursts were observed 
by RXTE and improves our estimates of the Eddington fluxes 
of the other sources listed. 

The Eddington fluxes for each source are calculated in the 
same way as in Appendix A of IWGP131 Where no uncer¬ 
tainty on a flux measurement is given in the literature, we have 
conservatively assumed it is 25% of the flux measurement. 
The addition of the supplemental bursts did no t sign ificantly 
change any of the Eddington fluxes listed in lWGP13l 

The number of PRE and non-PRE bursts listed in Table 
[H represent the bursts which made it into our final analysis; 
bursts discarded for source confusion or the inability to fit 
a persistent emission model are not counted. These are still 
suitable to estimate the Eddington flux of the source. For ex¬ 
ample, IGR 17473-2721 has three PRE bursts in the complete 
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• Eddington Flux 
■ Brightest non-PRE burst 
▲ Predicted F Edd lower limit 
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Fig. 9.— Lower limits of Fgdd for various neutron stars from their non-PRE bursts using the method described in 0 We have plotted the brightest non- 
PRE burst (squares) and our new prediction (triangles). For sources with known Fgdd these are shown (circles) to demonstrate that the new method does not 
overestimate the Eddington flux. For sources where Fe dd has not been directly measured, we have obtained a greatly better constraint for the Fgdd of Cir X-1 
and IGR 17480-2446, and a marginally better constraint for 1M 0836-425, than those given by their brightest non-PRE bursts. For points with no visible 1<t 
errorbars, these are smaller than the marker. 
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TABLE 3 

Supplemental PRE bursts 


Source 

Burst Time 
(MJD) 

Instrument 

Flux 

(10~ 9 erg s _1 cm -2 ) 

Reference 

GX 17+2 

44749.0 

Hakucho 

9.3 ± 1.0 

i 

GX 17+2 

44749.5 

Hakucho 

12.2 ± 1.1 

i 

GX 17+2 

45194.8 

Hakucho 

10.1 ± 1.1 

i 

GX 17+2 

45197.6 

Hakucho 

12.8 ± 1.0 

i 

4U 2129+12 

47454.7 

Ginga 

40.0 

2 

SLX 1735-269 

52897.7 

JEM-X 

57.6 ± 9.6 

3 

SLX 1735-269 

53108.0 

IBIS/INTEGRAL 

64.0 

4 

SLX 1735-269 

56267.2 

MAXI/GSC 

48.0 ± 3.0 

5 

GX 3+1 

53248.8 

JEM-X 

53.4 

6 

4U 1722-30 

42666.0 

OSO-8 

60.0 

7 

4U 1746-37 

46317.0 

EXOSAT 

10.0 ± 1.0 

8 

4U 1746-37 

46318.0 

EXOSAT 

9.0 ± 1.0 

8 

4U 1820-303 

42550.0 

SAS-3/RMC 

39.0 ± 2.0 

9 

4U 1820-303 

42940.0 

SAS-3/RMC 

53.0 ± 4.2 

9 

4U 1820-303 

42941.0 

SAS-3/RMC 

41.0 ±4.5 

9 

4U 1820-303 

42942.0 

SAS-3/RMC 

46.0 ± 3.7 

9 

4U 1820-303 

42942.0 

SAS-3/RMC 

43.0 ±3.9 

9 

4U 1820-303 

42944.0 

SAS-3/RMC 

42.0 ± 2.9 

9 

4U 1820-303 

46296.0 

EXOSAT 

52.8 ± 1.9 

10 


REFERENCES. — 1 . iTawara et alj 419841). 2. Ivan P aradiis et al. 1 1990), 3. Molkov eL al]j2005l). 
4. iNego ro et al. 1 2012), 5. Sguera et al. 1 2007). 6. Chenevez et al. 1 2006), 7. Swank et alj 4 19771) , 8. 
Sztajno et al. i 1987), 9. Vacca et al. 1 1986), 10. Damen et al. 1 1990) 
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RXTE catalog and therefore has a known /*Edd- but all three of 
these bursts are source confused so they are not investigated 
in this work or listed in Table [j] 
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